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Using a hydrodynamic model, we predict the transverse 
momentum dependence of the spectra and the elliptic flow 
for different hadrons in Au+Au collisions at y^s = 130 A GeV. 
The dependence of the differential and p t -integrated elliptic 
flow on the hadron mass, equation of state and freeze-out 
temperature is studied both numerically and analytically. 

PACS numbers: 25.75-q, 25.75.Ld 
Keywords: Relativistic heavy-ion collisions; Elliptic 
flow; Hydrodynamic model 

1. Introduction. - One of the first observables measured 
at the Relativistic Heavy Ion Collider (RHIC) was the 
so-called elliptic flow jij. It describes the azimuthal mo- 
mentum space anisotropy of particle emission from non- 
central heavy-ion collisions in the plane transverse to the 
beam direction. Elliptic flow is characterized by the sec- 
ond harmonic coefficient V2{y, pt) of an azimuthal Fourier 
decomposition of the momentum distribution m . We 
here discuss elliptic flow at midrapidity, y = 0; its pt- 
averaged value is denoted simply by i> 2 . 

Elliptic flow is a fundamental observable since it di- 
rectly reflects the rescattering among the produced parti- 
cles. Rescattering transfers the initial spatial anisotropy 
of the nuclear overlap region in the transverse plane to 
the observed momentum distribution. For a given initial 
spatial deformation, the largest elliptic flow coefficient is 
obtained in the hydrodynamic limit where rescattering is 
so intense that the matter in the reaction zone reaches 
a state of local thermal equilibrium. Since the spatial 
anisotropy is largest at the beginning of the evolution, el- 
liptic flow is especially sensitive to the early stages of sys- 
tem evolution |5|(|. A measurement of v 2 thus provides 
access to the fundamental thermalization time scale in 
the early stages of a relativistic heavy- ion collision . 

In a preceding Letter we showed |J that the first 
measurements at RHIC of the elliptic flow of charged 
particles can be satisfactorily described by a hydrody- 
namical model, with initial and freeze-out conditions ob- 
tained by a straightforward extrapolation from a simi- 
lar analysis of Pb+Pb collisions at the SPS @J|. The 
data deviate, however, from the hydrodynamical predic- 
tion at large impact parameters b^7 fm, where the reac- 
tion volumes become small, and at large transverse mo- 
menta pt^tl-5 GeV/c. At the lower SPS energies (y/s = 
17 A GeV) the available data 0|nj are less conclusive, 



but indicate discrepancies with the hydrodynamic ap- 
proach for all but nearly central collisions |],[l2| ■ Whereas 
in the deviations from hydrodynamic predictions 

were interpreted to signal incomplete thermalization, the 
authors of |L2| suggest (at least in connection with the 
SPS data) that they may be due to the transition from 
an early hydrodynamic to a late hadronic kinetic stage, 
modeled by a hadronic cascade (RQMD) which leads to 
an earlier saturation of both radial and elliptic flow. 

In the present paper we present hydrodynamic pre- 
dictions at y/s = 130 A GeV for the shape and impact 
parameter dependence of the single particle spectra and 
for the magnitude and pt-dependence of the elliptic flow 
for a variety of hadron species. Such data should soon 
become available. They will not only provide evidence 
for the degree of thermalization of other hadronic species 
and thus provide valuable input for the discussion of the 
two different interpretations just mentioned |],|l2|], but 
also help to select between the different combinations of 
initial and freeze-out conditions and equations of state 
studied in g|. 

2. Hydrodynamic Results.- The hydrodynamic code, 
its initialization and the calculation of final state spec- 
tra and elliptic flow coefficients are described in [^|,[|. 
We here use the parameter sets given in the last three 
columns of Table I in and label the results accord- 
ingly by EOS Q(120), EOS Q(140), and EOS H(140), 
respectively. EOS H is a hadron resonance gas equation 
of state with sound velocity c s « V0.15; EOS Q includes 
a first order phase transition to an ideal quark-gluon gas 
(c s = y/TJl) at T c = 164 MeV with a latent heat of 
1.15 GeV/fm 3 |jl3"|| . In the energy region studied here, 
EOS Q is effectively softer than EOS H, since the ex- 
pansion is largely controlled by the soft phase-transition 
region ||. Freeze-out is implemented along a surface of 
constant energy density; the numbers in brackets denote 
the approximate freeze-out temperature in MeV on this 
surface. 

Radial flow and single-particle spectra: The transverse 
momentum spectra for negative hadrons, neutral pions 
and protons are shown in Figure [l]. All decay prod- 
ucts from strong or electromagnetic decays of unstable 
resonances up to 1.4 GeV mass are included. These pt- 
spectra are considerably flatter than the corresponding 
SPS spectra; the thin solid lines show our earlier fits 
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to the data from Pb+Pb collisions at y/s = 17 ^4 GeV 
Ji"4j . This reflects the stronger radial flow at RHIC which 
hydrodynamics predicts as a result of the higher initial 
energy density |^ . The crosses in Figure [j] show that the 
hadronic cascade code UrQMD |H| gives much steeper 
spectra [jl6); in fact, the UrQMD spectra at RHIC ener- 
gies are slightly steeper than the measured spectra at SPS 
energies. Apparently the scattering mechanisms built 
into UrQMD are not efficient enough to build a suffi- 
cient amount of radial flow. On the other hand, the flat- 
ter slopes of hydrodynamic pion spectra are consistent 
with first preliminary data on negative hadron spectra at 
midrapidity jl7) . The final data are expected to be accu- 
rate enough to distinguish between the three parameter 
sets shown in Figure [l] (or exclude all of them) . 
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FIG. 1. The p t spectra of neutral pions (upper left), nega- 
tive hadrons (upper right), and protons (lower left) for the 6% 
most central Au+Au collisions at ^J~s — 130 A GeV, for differ- 
ent equations of state and freeze-out temperatures (see text). 
The crosses show UrQMD results for neutral pions and pro- 
tons at ^/s = 200 A GeV @ for comparison. Also shown are 
the corresponding hydrodynamic spectra at ^J~s = 17 A GeV 
for EOS Q(120) (thin solid lines). 

Figure ^| shows mt-spectra for pions, kaons, protons 
and f2 hyperons at different collision centralities. For 
peripheral collisions the spectra get steeper, due to a de- 
crease of the average radial flow velocity (v±) as a result 
of a lower initial energy density and a shorter lifetime of 
the reaction zone. The comparison between theory and 
experiment will be particularly interesting in the region 
of large impact parameters where the measured elliptic 
flow Q lags behind the hydrodynamical results |)|. As 
discussed in the Introduction, the single-particle spectra 
may help to understand whether the origin of this dis- 
crepancy is incomplete thermalization at an early or late 
stage of the evolution. 
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FIG. 2. The mt-spectra of negative pions (upper left), 
kaons (upper right), protons (lower left) and Q baryons (lower 
right) for Au+Au collisions at yfs = 130 A GeV, for colli- 
sion centralities (top to bottom) b < 5.4, 5.4 < b < 9.9 and 
9.9 < b < 13.5 fm. For clarity the spectra for different cen- 
trality bins are separated by factors of 10. The calculations 
were done with EOS Q. The fi distribution is also shown for 
Tf = 164 MeV to simulate decoupling at the phase transition. 
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Table 1. Elliptic flow coefficient V2 at midrapidity for various 
hadron species from minimum bias Au+Au collisions at ^J~s — 
130 A GeV. The pt-average was taken over the full spectrum. 



Elliptic flow. The impact parameter and pt depen- 
dence of the elliptic flow coefficient vi for charged par- 
ticles from Au+Au collisions at y/s — 130 A GeV was 
discussed in B. Table 1 shows predictions for the pt- 
integrated elliptic flow of a variety of identified hadrons 
in minimum bias collisions, for the same three parame- 
ter sets as studied in that paper. For the <p meson and 
the fl hyperon we include two values, one for simulta- 
neous freeze-out with all other hadrons, the other for 
freeze-out directly after hadronization at T c = 164 MeV. 
The latter option accounts for the expectation that the 
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heavy and weakly coupled fl is not likely to participate 
after hadronization in any modifications of the previ- 
ously established flow pattern JOl. Similar arguments 
hold for the (f> meson. Even though at RHIC energies 
a large fraction of the flow is already established be- 
fore hadronization ||, the additional flow generated af- 
terwards by hadronic rescattering is seen to affect the 
f2 quite strongly, for both w 2 (Table 1) and the single- 
particle slope (Figure ||). 

For identified pions the dependence of «2 on EOS and 
Tf is similar as for the sum of all charged hadrons (see 
) : Lower freeze-out temperatures and harder EOS lead 
to flatter single particle spectra and thereby to larger p t - 
integrated elliptic flow. For identified hadrons we see 
that, as their mass increases, there is stronger sensitivity 
to the EoS than that to Tf (see also Figure || below) . 
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FIG. 3. pt-differential elliptic flow at midrapidity for vari- 
ous hadrons from minimum bias Au+Au collisions at ^/s — 
130AGeV for EOS Q(120). 

Figure [| shows the differential momentum anisotropy 
V2(pt) for different hadron species for EOS Q and Tf ss 
120 MeV. At a given value of pt, the elliptic flow is seen to 
decrease with increasing particle mass. This is a conse- 
quence of rest-mass-dependent radial flow effects on the 
shape of the single-particle pt-spectrum, as will be ana- 
lytically discussed in the following section. 

The smaller differential anisotropy at fixed p t does not 
contradict the results in Table 1 which generically give 
larger p t -averaged elliptic flow for heavier particles. This 
is a consequence of the fact that radial flow leads to a 
flattening of the pt-spectra of heavier particles p9| , p0[ , 
whose pt-averaged vi thus receives more weight from the 
high-pi region where V2(pt) is larger. Whether this larger 
spectral weight for high pt wins over the reduction of 
vi at fixed p t depends on the details of the expansion 
dynamics. 

The effect of the EOS and the freeze-out temperature 
on the differential elliptic flow of pions and protons is 



demonstrated in Figure [J. The EOS affects V2(pt) for 
all hadrons in the same way: the stiffer EOS H leads to 
larger w 2 at low p t and to smaller v 2 at high p t than the 
softer EOS Q. The effect of the freeze-out temperature 
on V2(pt) is more delicate: for pions the effect is small, 
and for EOS Q an increase in the freeze-out tempera- 
ture decreases both the p*-averaged and the differential 
elliptic flow. The heavier protons behave oppositely: the 
differential anisotropy increases with increasing freeze- 
out temperature. The origin of this behaviour will be 
studied in the following section. Clearly, the different 
Tf-dependence of ^(p*) of different particles can be used 
to constrain the freeze-out temperature independently of 
the EOS. 
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FIG. 4. The effect of the EOS and the freeze-out tem- 
perature on the elliptic flow of midrapidity pions (left) and 
protons (right) from minimum bias Au+Au collisions at 
■sTs = 130 A GeV. 

In view of the observed deviations from hydrodynamic 
behaviour of charged particle elliptic flow at large impact 
parameters Q it will be interesting to study the central- 
ity dependence of i>2 separately for a variety of hadron 
species. Corresponding hydrodynamic predictions are 
shown in Figure ^ (again including the option that the f2 
freezes out early at Tf = T c ). Particle-specific deviations 
from these predictions should provide valuable insights 
into the thermalization and freeze-out mechanisms. 

3. Analytical results.- In the remainder we try to un- 
derstand the hydrodynamic behaviour of V2(pt) and its 
dependence on the hadron mass and freeze-out tempera- 
ture, using a simple analytical model. Before going into 
the technical details we give a simple intuitive argument 
why, at small p t , the elliptic flow of heavier particles is 
smaller than for lighter ones. It is well-known that radial 
flow shifts the p t -distributions to larger values of pt, and 
that for nonrelativistic pt < m this effect increases with 
the particle mass m and the radial flow velocity (v±). In 
the extreme case of a thin shell expanding at high veloc- 
ity, the spectrum actually develops a relative minimum 
at pt = and a peak at nonzero p t ( "blast wave peak" 
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Jl9[), and with increasing mass and (v±) the peak shifts 
to larger p t . Relative to the case without radial flow, the 
spectrum is thus depleted at small p t , and the depletion 
as well as the pt range over which it occurs increase with 
m and (v±). 
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FIG. 5. The p t - differential elliptic flow of negative hadrons 
(upper left), pions (upper right), protons (lower left) 
and omega baryons (lower right) for Au+Au collisions at 
y/s = 130 A GeV for collision centralities (top to bottom) 
9.9 < b < 13.5, 5.4 < b < 9.9 and b < 5.4 fm. The cal- 
culations were done for EOS Q. 

If, as in the case of fully developed elliptic flow, the 
radial velocity is larger in x than in y direction, \v x \ > 
\v y \, the same is true for this relative depletion effect. It 
counteracts the excess of particles with p t in a;-direction 
over those with pt in y-direction which is the origin of 
t>2 > 0. Thus it reduces V2 at small pt- This reduction 
and the pt-range over which it occurs both increase with 
the particle mass m and the average radial flow (v±). 
A quick graphic sketch shows that in the extreme case, 
where the single-particle spectrum develops a "blast wave 
peak" ]l9| , w 2 even turns negative at low p t , due to the 
shift of the peak to larger p t in x than in y direction if 

\V X \ > \Vy\. 

When looking for a simple model which allows to show 
this analytically, a generalization of the "blast wave" 
model fll9fl comes to mind, in which thermalized mat- 
ter of temperature Tf , approximated by a boosted Boltz- 
mann distribution, freezes out on a thin cylindrical shell 
along which the radial flow shows an azimuthal modula- 
tion with \v x \ > \v y \. Assuming boost-invariant longitu- 
dinal expansion and freeze-out at constant proper time, 



the Cooper-Frye freeze-out spectrum can be calculated 
by trivially generalizing the derivation given in pl| , p2| . 
In the Boltzmann approximation^ one finds, up to irrel- 
evant constants, 



dN 



dy dm\ < 
where <£ 



2tt 



# s tfi (&(&)) e Qt(W cos ^-"« , (I) 



are the azimuthal angles in coordinate 



and momentum space, and the arguments at(4> s ) = 
(p t /T)smh(p(<f> s )), /3 t (^)=(m t /T) cosh(pf> s )) §§ now 
depend on a ^-dependent radial flow rapidity p{4> s )- 
The elliptic flow coefficient is obtained by taking the az- 
imuthal average over cos(2</> p ) with this spectrum, V2 — 
(cos(2cf>p)} . The </> p -integral can be done analytically: 

f^defts cos(2^)/ 2 (a t (^))ifi(/3 t (^)) 

V2\Pt) = — ^ ■ (2) 

r o 27r d0 s / o (at(^))^i(A(^)) 

Making the Ansatz p = po + p a cos(20 s ) we checked that, 
for a suitable choice of the average radial flow rapidity 
Po and its azimuthal modulation amplitude p a , Eq. (Q) 
reproduces all relevant features of the full hydrodynamic 
results almost quantitatively. The (/^-independent ver- 
sion of Eq. (|l|) |^(| has been used very successfully by 
many groups to fit single-particle spectra from central 
heavy-ion collisions in order to extract the average tem- 
perature and radial flow velocity at freeze-out. Similarly, 
Eqs. (|l|) and (||) can be used to fit the spectra and elliptic 
flow from non-central collisions in order to extract also 
the average azimuthal flow modulation. 



FIG. 6. Simple source of four fireballs. 

Unfortunately, the remaining angular integral in (|^) 
cannot be done analytically. To reach a fully analytical 
understanding of many of the features observed in the 
previous section one can exploit an even simpler model 
which still captures the main effects qualitatively, but no 
longer quantitatively. It consists of four non-expanding 
fireballs of equal volume, freezing out instantaneously 



1 For Bose-Einstein or Fermi-Dirac distributions, one simply 
substitutes Eq. (jjj) by a sum of terms with T replaced by T/n, 
n = 1, 2, . . ., and weighted with (±1)" [|l). 
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and simultaneously in the laboratory frame, whose cen- 
ters move in the transverse plane without longitudinal 
velocity component (see Fig. |J). What matters for the 
observed particle flow pattern are the velocities of the 
four fireballs, as indicated in the figure, but not their 
positions nor the size of their common laboratory frame 
volume. 

For this model one finds, again in the Boltzmann ap- 
proximation, the simple expression 

r (»)_ e #( 7 ,- 7B ) j2 («) 

V2{y ' Pt) = M^pj +e* .)J (™)' ( ) 

In this equation the particle mass enters only in the term 
e —hx-i v )^ an d it j s eaS y to see that if all other variables 
are held fixed, vi decreases with increasing mass. Fig- 
ure shows that our schematic source also reproduces the 
feature of negative V2 for protons at small pt which was 
anticipated above from the "blast wave" model. (This 
feature is also preserved by the cylindrical shell model, 
Eq. (||).) The "blast wave peak" in the single-particle 
spectrum is known to disappear when the constant ex- 
pansion velocity is replaced by a realistic radial velocity 
distribution ; we therefore then also expect the dip of 
V2 to negative values to be washed out, in agreement with 
Figs. H and EJ. On the other hand, the UrQMD calcula- 
tions reported in [ p3| do show negative V2 for midrapidity 
protons at small pt ; it would be interesting to further in- 
vestigate the origin of this effect in that model. It may 
be worth pointing out that a similar interplay between 
radial flow, random thermal motion and a directed flow 
anisotropy (first instead of second harmonic modulation) 
has been shown p3] to be responsible for negative di- 
rected flow coefficients v\ at low p t [Eol . 
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FIG. 7. Transverse momentum dependence of elliptic flow 
for midrapidity pions and protons from the schematic source 
in Figure |, for T = 140 MeV, v x = 0.6, and v y = 0.5. 

Our schematic source also permits us to understand 
the approximately linear p t -dependence of v% in the in- 
termediate pt region. Expanding the Bessel functions in 



(||) for large arguments and keeping only the leading term 
gives for midrapidity (y=0) particles 



t>2 (pt) ~ tanh 



1 / np t - Xm t 

2 V T 



(4) 



where k 



lxv x - jyv y , \ = j x — J y and fi = ln(J^). 

Since k > A for \v x \ > \v y \, we get for p t ^> m the simple 
form V2(pt) ~ tanh(^p t ) with ^ = (n — X)/2T. For the 
example in Fig. [7], £ « 5 GeV, and V2 begins to turn 
over at pt > 2.5 GeV, saturating at D2 = 1 as pt — > 00. 
This consideration also shows that the scale at which V2 
changes from its required quadratic pt-dependence at low 
Pt p^j (with a positive or negative coefficient) to a linear 
behaviour at intermediate p t is given by the particle rest 
mass m, multiplied by a prefactor of order 1 which is 
governed by the average radial flow velocity (v±). 

We finally discuss the dependence of V2 on the freeze- 
out temperature. Taking the temperature derivative of 
Eq. (Eh we get (again at y=0) 



dv 2 , s _ Am t 
dT ^ Pt ' ~ 2T 2 



KPt 



1 



cosh 



T 



(5) 



With k > A one sees that at fixed pt the sign of the deriva- 
tive depends on the particle mass and that, for increasing 
Tf , v-iijPt) may thus decrease for pions while increasing for 
protons, as seen in Fig. Ej^ 

4- Summary.- We have presented a variety of pre- 
dictions for the elliptic flow and single-particle spectra 
for different hadron species produced in Au+Au colli- 
sions at y/s = 130^4 GeV, using a relativistic hydrody- 
namic model. We have studied the sensitivities to the 
equation of state and freeze-out temperature and showed 
that these can be used to further constrain the model pa- 
rameters and test the approach on a quantitative level. 
A simple expression for fitting spectra and elliptic flow 
data in order to extract the average radial flow and flow 
anisotropy has been given. Crucial features of the p t - 
dependence of the elliptic flow have been elucidated with 
a simple schematic model. Testing the predicted pt- 
dependence of i>2 for many different hadron species will 
clarify the validity of the picture of a thcrmalizcd expand- 
ing source with a common flow velocity for all hadrons 
at RHIC energies. 
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2 Of course, in a full hydrodynamic calculation Tf cannot be 
varied independently of the flow velocities, and the latter are 
not fixed, but vary over the freeze-out surface. Eq. (|^) can 
only give a qualitative impression. 
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